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The problem of Kramers' turnover is a central issue of dynamical theory of reaction rate. Since 
its classical solution in the Markovian limit in mid-eighties by Melnikov and Meshkov , the problem 
has been addressed by a number of groups in the last decade both in classical non-Markovian and 
quantum mechanical context. Based on a coherent state representation of noise operators and a 
positive definite Wigner canonical thermal distribution function we have recently developed a c- 
number quantum Langevin equation [Barik et al, J. Chem. Phys. 119, 680 (2003); Banerjee et 
al, Phys. Rev. E 65, 021109 (2002)]. We implement this scheme within PoUak's weU known 
normal mode description to calculate the quantum transmission coefficient over an arbitrary range 
of friction, noise correlation and temperature. The theory generalizes the quantum correction to 
Grote-Hynes factor in the rate expression down to vacuum limit which reduces to well known high 
temperature quantum correction, i.e., the Wolynes term for quantum transmission and reflection for 
the barrier in the appropriate limit and also considers the quantum corrections due to nonlinearity of 
the system potential order by order which contributes to energy loss and dispersion due to coupling 
between unstable and stable normal modes near the barrier top and is valid for both above and 
below the activated tunneling regime. Our results have been compared with those obtained earlier 
for a model potential and found to be good agreement. 

I. INTRODUCTION 

Kramers' 0] diffusion model of chemical reactions proposed in 1940 forms the dynamical basis of modern rate 
theory of activated processes. The seminal and essential content of this nonequilibrium formulation is the inclusion 
of dependence of rate constant of a reaction on viscosity or friction of the reaction medium. Based on the classical 
theory of Brownian motion in phase space Kramers derived the expressions for nonequilibrium steady state distribution 
functions to work out the rate coefRcients in the two different limiting situations and showed that the rate varies linearly 
in the weak dissipation regime and inversely in the high dissipation regime with friction. That is, in between the 
energy diffusion and spatial diffusion limited regimes the rate constant as a function of friction exhibits a bell-shaped 
curve - known as Kramers' turnover. With the advent of ultrafast lasers and timc-rcsolved detection techniques 
since late seventies, experimental confirmation of Kramers' theory provided a new impetus for further development 
in chemical dynamics and condensed matter physics H, H, IE 111 - 

In spite of this spectacular growth one point however, that remained illusive for several decades was the absence 
of an analytic rate formulae which correctly interpolates the two limits. A significant advancement was made in 
1986 when Melnikov and Meshkov (71, |8| proposed a theory of Kramers turnover for Markovian friction. The key 
point of this analysis is the calculation of an average energy loss of the system due to its weak coupling with the 
bath modes near the top of the potential barrier. The approach was further extended to non-Markovian domain by 
PoUak, Grabert and Hanggi 9] within the framework of a normal mode description of the system plus reservoir, where 
the implementation of an ad hoc function in the expression for the spatial diffusion limited rate could be avoided. 
Subsequent to this formulation of classical theory, the problem of quantum Kramers' turnover was addressed by Rips 
and PoUak ^3 ■ The degree of accuracy of classical and quantum theories has been tested by numerical simulation of 
reaction rate with model potentials. We refer to Refs. 11-16 for further details. 

The classical theory of Kramers' turnover uses a Langevin description that governs the dynamics of the system 
by an infinite number of harmonic oscillators coupled linearly to the system degree of freedom. Very recently based 
on a coherent state representation of noise operator and a positive definite Wigner canonical thermal distribution of 
harmonic oscillators of the bath we have proposed a c- number quantum Langevin equation [l^lT9ll20ll2lIl2^ l2^ 
in the context of quantum rate theory and stochastic processes in terms of phase space function formalism |24L l25l | . It 
would seem that one should be able to analyze a quantum turnover for arbitrary noise correlation and temperature 
by using such a c-number formulation. This is the main purpose of this paper. In what follows we make use of a 
c-number Hamiltonian in the normal mode procedure and take care of the nonequilibrium dynamics at the barrier 
top by calculating the average quantum energy loss of the unstable mode due to its weak coupling with the stable 
bath modes and cquilibriation in the well by a Wigner thermal distribution. Besides being an approach based on 
canonical quantization the theory specifically addresses the following three points: First, we quantize the system 
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mode to make the theory apphcable beyond the activated tunnehng regime down to absolute zero. Second, although 
harmonic approximation is a good description at the barrier top, actual shape and nonlinearity of the system potential 
is important at low temperatures; but systematic corrections to harmonic approximation have hardly been envisaged. 
The present formulation takes care of quantum corrections due to nonlinearity of the system degrees of freedom order 
by order. Again even in the lowest order this quantum correction can be, in principle, a contribution in the weak 
coupling between the stable and unstable modes for calculation of energy loss and dispersion as shown in this paper. 
Third, we take into account of the quantum correction to Grote-Hynes dynamical factor in the rate expression down 
to vaccum limit. This reduces to well-known high temperature quantum correction or the Wolynes |26l | term in the 
appropriate limit. The approach thus generalizes the Wolynes term in the deep tunneling regime. 

The outlay of the paper is as follows. In the following section II we discuss our c- number quantum Langevin equation 
which allows us to realize a c-number Hamiltonian analyzed by PoUak's normal mode analysis. The equilibrium theory 
in terms of Wigner distribution function to formulate a quantum counterpart of multidimensional transition state 
theory (TST) has been presented in Sec. III. Sec. IV is devoted to nonequihbrium dynamics of the unstable mode and 
its average energy loss over round trip time in the well near the barrier top. Since this energy is sensitive to quantum 
contribution due to nonlinearity of the potential we calculate the quantum correction to energy loss of the unstable 
mode explicitly in Sec.V. An explicit example with cubic potential has been worked out in Sec. VI to compare with 
the known results. The paper is concluded in Sec. VII. 

II. C-NUMBER QUANTUM LANGEVIN EQUATION 

We consider a particle of unit mass coupled to a medium comprised of a set of harmonic oscillators with frequency 
LUi. This is described by the following Hamiltonian: 

i=l ^ ■> 

Here q and p are co-ordinate and momentum operators of the particle and the set {xi,pi} is the set of co-ordinate 
and momentum operators for the reservoir oscillators coupled linearly to the system through their coupling coefficients 
Ci. The potential V{q) is due to the external force field for the Brownian particle. The co-ordinate and momentum 
operators follow the usual commutation relations [(1,01 — ih and — ihSij. 

Eliminating the reservoir degrees of freedom in the usual way we obtain the operator Langevin equation for the 
particle, 

m + f dt'jit - t')'q{t') + V'{q) = F{t) , (2.2) 

JO 

where the noise operator F{t) and the memory kernel 7(t) are given by 



3 

and 

N ^2 

Here masses have been assumed to be unity. 

The Ea. (|2.2|l is the well known exact quantized operator Langevin equation for which the noise properties of F{t) 
can be derived by using a suitable initial canonical distribution of the bath co-ordinate and momentum operators at 
i = as follows; 



1 ^ 
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{F{t))Qs = (2.5) 
\{mF{t'))Qs + {F{t')F{t))Qs} (2.6) 



2 

N 



where {...)qs refers to quantum statistical average on bath degrees of freedom and is defined as 

{d)Qs = T'^Q^^p(-^^'^ath/fcBr) ^2.8) 

Trexp(-i/bath/fcs7') 

for any operator d{{{uj] /cj)xj - q},{pj}) where -ffbath(I]jIi(Pi /2 + l/2{uJiXi - ^AY)) at i = 0. By Trace we 
mean the usual quantum statistical average. Ea. H2.7|l is the fluctuation-dissipation relation with the noise operators 
ordered appropriately in the quantum mechanical sense. 

To construct a c- number Langevin equation we proceed from Eq. (|2.2|) . We first carry out the quantum mechanical 
average of Eq. H2.2|) 



(m + f dt'^{t - t'mt')) + {V'{q)) = (Fit)) (2.9) 
Jo 

where the quantum mechanical average (. . .) is taken over the initial product separable quantum states of the 
particle and the bath oscillators at t = 0, |(/)){|Q!i)|a2) . . . |aAr)}. Here |0) denotes any arbitrary initial state of 
the particle and \ai) corresponds to the initial coherent state of the i-th bath oscillator. \ai) is given by \ai) — 
exp(— |aip/2) ^^^p(a"'/VniT)|rii), ai being expressed in terms of the mean values of the shifted co-ordinate and 

momentum of the i-th osciUator, {{ujf /ci){xi{0)) ~ (g(0))} — {Vh/2L0i){ai + a*) and (pi(0)) = iy^huji/2{a* — at), 
respectively. It is important to note that {F{t)) of Eg. 12.91) is a classical-like noise term which, in general, is a non-zero 
number because of the quantum mechanical averaging and is given by {{F{t)) = f{t)); 



m = E 



— (%(0)) - (9(0)) \^cosu;jt+^{p,{0))smcj,t 



(2.10) 



It is convenient to rewrite the c- number equation (|2.9|l as follows; 



m + / dt'j{t ~ t'mt')} + {v'm ^ nt) (2.11) 







To realize f{t) as an effective c-number noise we now introduce the ansatz that the momenta {pj{0)) and the shifted 
co-ordinates {{uj'j /cj){xj{0)) — (g(0))}, {pj} of the bath oscillators are distributed according to a canonical distribution 
of Gaussian form as 



Fj = Nexp< 



3 

2hwj (fijiujj) + i) 



(2.12) 



so that for any function of the quantum mechanical mean values Oj{(pj(0)), 
{{uj'j /cj){xj{0)) — (g(0)))} the statistical average (. . .)s is 



{Oj}s - I O, P, d{p,{0)) d{(c.2/c,)(£,(0)) - (q(0)>} 



(2.13) 
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Here fij indicates the average thermal photon number of the j-th oscillator at temperature T and nj{LUj) = 
l/[exp {hwj/kBT) — 1] and N is the normalization constant. 

The distribution (j2.12|l and the definition of statistical average H2.13II imply that f{t) must satisfy 



(/(i)>5 = (2.14) 

and 

{f{t)f{t'))s = \Y.^^^ ('°*^2Sr) '^°^^^(*-^') (2.15) 



That is, c- number noise f{t) is such that it is zero-centered and satisfies the standard fluctuation-dissipation relation 
(FDR) as expressed in Eq.j^TJ. It is important to emphasize that the ansatz H2.12|l is a canonical Wigner distribution 
for a shifted harmonic oscillator which remains always a positive definite function. A special advantage of using 
this distribution is that it remains valid as pure state non-singular distribution function at T = 0. Furthermore, 
this procedure allows us to bypass the operator ordering prescription of Ea. H2.7|l for deriving the noise properties of 
the bath in terms of fluctuation-dissipation relation and to identify f{t) as a classical looking noise with quantum 
mechanical content. 

We now return to Ea. l|2.11() to add the force term V'{{q)) on both sides of Ea. H2.11|l and rearrange it to obtain 

q = p (2.16) 

P = - I dt'j{t-t')p{t')-V'{q) + f{t) + Q{t) (2.17) 
Ja 

where we put {q(t)) — q{t) and {p{t) — p{t) for notational convenience and 



Q{t) = V'm-{V'{q)) (2.18) 

represents the quantum correction due to the system degrees of freedom. Eg. (|2. 17(1 offers a simple interpretation. 
This implies that the quantum Langevin equation is governed by a c- number quantum noise /(<) originati ng fr om the 
heat bath characterized by the properties H2.14|l and 12.15|l and a quantum fluctuation term Q{t) [islll^ I20ll2lll23 . l2^ 

characteristic of the non-linearity of the potential 29] . 

Referring to the quantum nature of the system in the Heisenberg picture, one may writes. 



q{t) = q + Sq (2.19) 
Pit) = p + 5p (2.20) 

where {q){= q) and (p)(= p) are the quantum-mechanical averages and Sq, dp are the operators. By construction 
{6q) and (Sp) are zero and [(5q, ^p] = ih. Using Eos. (|2. 19(1 and (|2.20() in {V'{q)) and a Taylor series expansion around 
{q) it is possible to express Q{t) as 

Q{t) = -Y,-V^-+^Hq){Sq"{t)) (2.21) 
n>2 

Here V^'^\q) is the n-th derivative of the potential V(q) ■ For example, the second order Q{t) is given by Q{t) = 
-^V"'{q){6f). The calculation of Q{t) M, M, 2^^^ M therefore rests on quantum correction terms, 
((5q"(<)) which are determined by solving the quantum correction equations as discussed in the Sec.V. 

The c-number Hamiltonian corresponding to Langevin equation H2.16I I2.17|l is given by 



H 



EL 
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(2.22) 
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Note that the above Hamihonian is different from our starting Hamiltonian operator H2.1|l because of the c-number 
nature of H2.22I) . {xi,Pi} are the quantum mean value of the co-ordinate and the momentum operators of the bath 
oscillators. 

The spectral density function is defined as 



N 2 

JM = 5 ^ ^ ^(c^ - c^O (2.23) 

2—1 

We now assume that at g = 0, the potential V{q) has a barrier with height such that a harmonic approximation 
around q = leads to 



Viq) = Vi-^Lulq^ + V2{q) (2.24) 

where ujf = V"{q) \q=o, refers to the second derivative of the potential V{q). ojb is the frequency at the barrier top 
and V2{q) is the non-Hnear part of the classical potential and is given by V2 = J2n>3 ;7! ^ap^ l9=o <?"• With Eg. p. 24(1 
the quantum correction part in the Hamiltonian Eq. 12.22(1 becomes 



E ^y^''Hq){Sr) - -fB2 + VM (2.25) 

ri>2 

where Bn = V3 — X]n>2 %T^~Bq^ ■ Note that we have introduced an approximation by putting a bar over 

quantum dispersion {Sq"{t)) to indicate its time average since we will be concerned here with the energy loss of the 
system mode averaged over one round trip time, i.e., the time required to traverse from one turning point of the 
potential well to another and back. We will discuss this averaging in greater detail in Sec.V. Putting ((2.24(1 and 1(2.25(1 
in the Hamiltonian 1(2.22(1 we obtain 

H = Ho + VN{q) (2.26) 
where we have decomposed the Hamiltonian in the harmonic part Hq and the anharmonic part V^iq) as 



Ho = 



Pi 
2 



E 



2 2 

-.^bq 



UOiXi 



-qf 



(2.27) 



and 



VN{q) = V2{q) + VM (2.28) 
and = V^-^Lol 

V2{q) and Vz{q) are therefore classical and quantum anharmonic contributions to total anharmonic part of the 
Hamiltonian. The separability of the c-number Hamiltonian in the quadratic and nonlinear parts allows us to make 
a normal mode transformation to convert the quadratic Hamiltonian into a diagonal form. The method of normal 
mode analysis has been used extensively by PoUak and co-workers in classical and quantum theories of activation, 
tunneling and dephasing. For details we refer to [Ref. 9,10,27-30]. 

Following PoUak, we diagonalizc the force constant matrix T of the Hamiltonian ((2.27(1 with the matrix U 



UT = X^U 

where U provides the transformation from old co-ordinates to the normal co-ordinates 



(2.29) 
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/ P \ 

yi 

2/2 



u 



I 1 \ 

Xl 
X2 



(2.30) 



The c-number Hamiltonian of the unstable normal co-ordinate is given by 



1 1 ^1 



(2.31) 



The eigenvalues Xf and are expressible in terms of the coupling constant of the system and the bath implicitly 
as follows: 



N 



N 



1+E 



i = 1,2. ..TV 



(2.32) 



(2.33) 



where H2.32|l and l|2.33(l correspond to normal mode frequencies of the unstable mode and the i-th bath oscillator 
respectively. 

The transformation H2.30|l implies 



N 

g = Moo P + ^ UjO Vj 



(2.34) 



and it has been shown |^ that uqo and Ujo may be expressed as 



^00 



N 



ftK + Ag)^ 



(2.35) 



and 



Making use of the spectral density function 
be written in the continuum limit as 



N 



and Laplace transformation of ^(t) Eq. 



(2.36) 

and Ea. H2.35|) may 



l + jiXb)/Xt 



(2.37) 



and 



2 

1 + - 

TT 



duj 



(A? 



(2.38) 



7 



The two important identities in relation to orthogonal transformation matrices and the associated frequencies may 
be noted here for the dynamics at the barrier top and at the bottom of the well; 



N 



N 



(2.39) 



and 



N 



N 



(2.40) 



Here ujq and Ag are the frequencies of the system at the bottom of the well in the original co-ordinate and normal 
co-ordinate respectively. Similarly corresponds to the normal mode frequencies of the bath oscillators coresponding 
to a normal mode Hamiltonian at the bottom of the well, 



Tj/ "'^ -'2 I 1 \ 2 /2 

" " 2^ 2 °^ 



1 



4 = 1 



(2.41) 



Here p' and p' are coordinate and momentum of system mode respectively and y[ and y[ are coordinate and 
momentum of ith bath oscillator respectively at the bottom of the well in the normal coordinates. Aq and are given 
by 



Ag = c^oV 



N 



1 + E 



A^ 



N 



1 + E 



i = 1,2. ..TV 



(2.42) 



(2.43) 



The description of a Kramers' turnover of the rate constant from low to high friction limit requires the coupling 
between the normal modes. An important parameter which is relevant for describing the coupling perturbatively has 
been defined by 



V"oo / V"oo 



W / X 2 
"00/ 



which can be expressed further in terms of the frequency dependent friction as 



(2.44) 



1 

'=2 



Afc dXb 

Furthermore for future use we also define a parameter 



(2.45) 



N / X 2 

UiO 



"00 



whose Laplace transform is given by 0, l3^ 



(2.46) 



2 2 
^00 * 



+ s 7(s) -uj^ s2 - XI 



(2.47) 
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III. C-NUMBER QUANTUM VERSION OF MULTIDIMENSIONAL TST 

To begin with we consider the particle to be trapped in a weU described by a potential V{q) depicted schematically 
in Fig.l. In the normal mode description of {N + 1) oscillators according to the Hamiltonian (|2.31|l the bath modes 
and the system mode are uncoupled. Considering the unstable reaction co-ordinate to be thermalized according to 
the Wigner thermal canonical distribution of N uncoupled harmonic oscillators plus one inverted we have 



exp 



1 X2 „2 



hXo{no{Xo) + i) 



N 



1 X2„,2 



(3.1) 



z is the normalization constant. As usual this can be calculated using the distribution function inside the reactant 
well. For this it is necessary to consider the normal mode Hamiltonian at the bottom of the well expressed as Hq in 
the Eg. H2. 41(1 . The corresponding distribution in the well is 



Weq = Z ^ exp 



„/2 ^ 



^Ao(no(Ao) + \) 



N 



J|exp 



2»i ' 2^^i 



Mi(ni(A,) + i) 



(3.2) 



which can be normalized to obtain 



Ao 



N 



A,, 



27rriAo(no(Ao) + i) fJi 2TThAi{n,{A,) + i) 
The identity relation H2.4()|l can be used to transform 1(3. 3|l to the following form 



(3.3) 



N 

n 



27r;iAo(no(Ao) + i) fj^ 2TThAi{n,{A,) + i) 
Putting Eq. 1(3.4(1 in Eq. 1(3.1(1 we obtain after integration over the stable modes 



(3.4) 



2TrhXa{na{Xo) + ^) 



Xb 

— X exp 



1X2^2 



ftAo(no(Ao) + i) 



where 



The total energy of the unstable mode is 



1^ hXi{ni{Xi) + ^) 
l\ M.(n,(AO + i) 



e = Ip' + v,^ Ixy 



(3.5) 



(3.6) 



(3.7) 



The prime quantity for determination of rate constant is the distribution of energy of the unstable mode. Thus 
going over to an energy space so that the co-ordinate p, p are transformed to t, E, respectively with unit Jacobian, 
the equilibrium distribution function ((3.51) is given by, 



2TThXo{no{Xo) + i) LOb 



X exp 



E 



fiAo(rto(Ao) + 5) 



(3.8) 



The above distribution is valid for the energy of the unstable mode E > as well as E < . 
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A. Quantum multidimensional TST rate 

The rate of activated barrier crossing in terms of the equilibrium probabihty becomes 

f{E) dE (3.9) 

yt 

As the unstable mode remains uncoupled from the stable modes the former mode behaves deterministically and 
the recrossing does not occur in this case. 

Making use of the distribution H3.8|l in (|3.9|) we obtain the rate constant 



i^TST = 7, X exp 



/iAo(no(Ao) + i) 



(3.10) 



The above expression corresponds to a quantum multidimensional transition state rate constant. This is central 
result of this section. Apart from usual Kramers'-Grote-Hynes term Ab/wf, and ujo/2-k, the term arising out of classical 
transition state result, Eq.(3.10) contains two important factors. First, an exponential Arrhenius term where the usual 
thermal factor fcsT is replaced by /iAo(no('^o) + 5) includes quantum effects due to heat bath at very low temperature. 
In the high temperature limit it reduces to ksT and one recovers the usual Boltzmann factor. This term is essentially 
an offshoot of a description of thermal equilibrium by a canonical Wigner distribution of harmonic oscillators heat 
bath. Second term x can be identified as the quantum correction to Grote-Hynes factor or more precisely a vacuum 
corrected generalized Wolynes contribution for quantum transmission and reflection for the finite barrier. While usual 
Wolynes term takes care of the quantum effects at the higher temperature the factor x incorporates quantum effects 
at arbitrary low temperature. In what follows we show that the usual Wolynes term can easily be recovered for x iii 
the appropriate limit. 

B. Derivation of Wolynes factor from x 

We begin by noting that n{x) in x which is given by 

^ HXi{ni{Xi) + i) 



hA,{ni{A,) + i) 



(3.11) 



is the Bose distribution n{x) = {e^^/^"'^ — 1) 1. Neglecting the vacuum contribution 1/2 from the terms like 
hx{ri{x) + ^) and keeping only the leading order quantum contribution in the thermal limit we obtain 



1 \ hx / . , hx 



-1 



nx[nix) + - ..-Uinl.^ (3.12) 



Therefore x reduces to S (say) the Wolynes factor 



N X^ sinh 



2k bT 



-1 



From identities H2.39|l and (|2.40l) it follows 



=1 A J sinh ^If- 



TV N 



-1 



(3.13) 



2—1 i—1 
N N 

nA.=^?n-' (3-15) 
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respectively, and we have 

N 



i=l 

Making use of the relation (|3.16|l in ^'6.l'S\i we obtain 



„ _ Xpujb TT- sinh(Mj/2fcBr) (-^ ^7] 



Furthermore (Xquji,) / (woXt) can be rewritten as ^ (liAb/fcaT) "which may be approximated in the form 

■ h{h\o/2kBT} 



K/^o)i^iS^@f?|. Eq.(3.17) then reduces to 



„ ^ cjfc sinh(^Ao/2fcaT) sinh(^A,/2fcBT) 

" Wo sinh(ftAb/2fcBr) 11 sinh(;iAi/2fcBr) ' 

This is the wellknown Wolynes |lOl l26j | expression derived in early eighties as a higher temperature equilibrium 
quantum correction to Kramers'-Grote-Hynes dynamical factor to Kramers' rate. Both x and the Wolynes factor 
become unity in the classical limit. We conclude by noting that unlike Wolynes factor S, x is valid below cross-over 
temperature. 



IV. THE COUPLING BETWEEN THE NORMAL MODES, ENERGY LOSS AND THE QUANTUM 

RATE 



The rate formula H3.9|) is quite general and requires the knowledge of f{E). In the last section the quantum 
version of Kramers'-Grote-Hynes estimation of this rate rests on the replacement of /(E) by an equilibrium Wigner 
distribution function in phase space. A determination of f{E) from the dynamics of the energy diffusion process, 
however is based on the solution of the integral equation originally formulated by Melnikov Q and subsequently by 
others [i[l3. 



fiE) = / dE' P{E\E') f{E') (4.1) 
Jo 

which implies that f{E) can be determined by a conditional probability function P{E\E') dE that a system escaping 
the barrier region with energy E' of the unstable mode makes a round trip of the barrier with an energy between E 
and E + dE. Deep inside the well, the strong coupling between the stable and unstable modes brings f{E) close to its 
equilibrium value Weq, whereas a weak coupling prevails where the energy is close to the barrier energy. A detailed 
consideration of this weak coupling between the stable and unstable modes results in the conditional probability 
function P{E\E') which when made use in H4.1|l gives the rate constant beyond the multidimensional TST limit. 

To proceed further we construct the equation of motion for the normal modes corresponding to the Hamiltonian 
(lOTl) using ijT^ and 

N 

P-^IP -"00 V^("oo/0 + ^Mio yi) (4.2) 

1=1 

N 

yi + >Hyt = -UioV^iuoQp + ^UMyi) (4.3) 

i=l 

where denotes the derivative of Vn{x) with respect to x. An important point is noteworthy. In contrast to 
classical theories where Vn is purely an anharmonic classical contribution, V2, the present treatment incorporates 
quantum effects due to nonlinearity of the system potential V3 entangled with quantum dispersion terms as given in 
Eq. (|2.28|l . Since the coupling between the stable and unstable normal modes is very sensitive to the small variation of 
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energy we expect the quantum effect to contribute significantly to the energy loss mechanism and hence the calculation 
of depopulation factor. This is an important point of departure from the earlier treatment of Rips and PoUak who 
had considered the unstable mode to be classical and neglected this contribution. The applicability of the theory is 
thus not restricted to socalled activated tunneling region only. 

Defining a smallness parameter gi as gi — Uio / uoo and ^ = J29i may write the zero order equation of motion 
for the unstable mode as 

P - p = -uoo V^("oo P) (4.4) 

and for the stable modes as 

y^ +>^fy^=9^C (4.5) 

where C is a forcing function given by 

C{t) = ~uoo V;,{uoo p) (4.6) 

The parameter e <C 1 is for weak coupling. The forcing function Q{t) can be determined by the unperturbed motion 
of the unstable mode p. The general solution of Eq. (|4.5I) is given by 

2;.(i) = y.(0)cosA,i+M sinA,t+ f dt' '-^^^^^^^-^ ({t') (4.7) 

\ JO 

As noted earlier that at equilibrium the initial distribution of the stable normal modes yi{0) and yi{0) is a thermal 
canonical Wigner distribution so that 

{y^ms = {MO))s - {y^{0)y^ms = 
{Uo?)s = \Uy^io)')s = (mx.,) + i") (4.8) 



The total c-number energy of the stable bath modes during the traversal of p mode over a round trip time t 



^ 1 1 



i=l 



can be calculated form 14.7|l and may be put in the form 

E = E' - AE + SE (4.10) 

E' = J2i Ei{0) where Ei{Q) is the initial energy of the i-th bath oscillator. The gain of energy by the stable modes 
is equal to the loss of energy of the unstable mode. This is characterized by a systematic and a stochastic contribution 
to energy loss of the unstable mode AE and 6E, respectively. These are given by 

Ai? = i V gf /*' dt /*' dt' cos[XS - t')] (it) at') (4.11) 
^ Jo Jo 

^ ftp 

6E^"^g, I dt[ j),(0) cos(A,t) - y,(0) A, sin(A,i) ] ({t) (4.12) 

While AE represents a systematic energy loss due to coupling, 5E is a measure of instantaneous loss induced by 
Gaussian fluctuation around (E) = E' — AE such that 
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{6E) = and (SE^) = D 



(4.13) 



Making use of the relation ()4.8|1 . {SE^) can be calculated as 



N 



i=l 



(SE^) =^g} j dt I dt' hXi 



"'0 



ni{Xi 



cos[A,(i-i')] m at') 



(4.14) 



In the Markovian limit H4.14|l reduces to 



(SE^) = 2 hXb 



nb + - 



AE 



(4.15) 



The mean and the second moments of the energy fluctuations lead us to the conditional distribution function of 
the Gaussian form as 



P{E\E') 



V2^ 



exp 



{E-E' + AEf 
2D 



(4.16) 



The Markovian limit of the above expression can be obtained by replacing D by H4.15|l . 

We now proceed to calculate the rate explicitly. As is wellknown that as probability function f{E) tends to its 
equilibrium limit Weq the rate T reduces to the rate calculated by multidimensional transition state theory. If, in 
general, they differ by a depopulation factor as defined in 



J. ^0 Xb 
It TT ^~ ^^P 

ZTT UJb 



Vi 



/iAo(no(Ao) + ^) 

so that by comparing (|4.17|) with H3.9|) to include all energies we write 



(4.17) 



27r 1 Wf, 

It ^ T- exp 

^0 X A& 



_ftAo(no(Ao) + 5) 
Introducing one-sided Fourier transformation ansatz as 



f{E) dE 



(4.18) 



and 



1 

X Ai 



exp 



Vi 



f{E) exp 



fi.Ao(no(Ao) + 5) 
2{iX+l/2)AE.E' 



D 



dE 



(4.19) 



1 UJb 

- exp 
X Ab 



Vi 



f{E) exp 



hXo{no{Xo) + ^) 

2{iX+l/2)AE.E' 



D 



dE 



(4.20) 



it is possible to transform H4.1|l into Wiener-Hopf [71 Il2. ISSL l3q | equation with symmetric kernel which can be solved 
by standard technique first suggested by Melnikov and Meshkov in the rate theoretical context and subsequently used 
by others. The ultimate expression of rate can be worked out as 
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X exp - 



ftAo(no(Ao) + i). 



X exp 



1 



oo 




2AE'^ (A'^ + 1/4) 
D 



)] 



dX 



(4.21) 



— OO 



[A'2 + 1/4; 



This is the central result of the paper. The quantum rate is a product of several terms. The first term corresponds to 
classical well frequency factor according to transition state theory. The second is the classical Kramers'-Grote-Hynes 
factor for arbitrary frequency dependent friction. The third term x as given by l|3.6|l corresponds to equilibrium quan- 
tum correction to Kramers'-Grote-Hynes factor which reduces to high temperature quantum correction or Wolynes 
term in the appropriate limit as shown in Sec. III. The fourth term refers to the Wigner's canonical thermal distribution 
in the harmonic well and reduces to the usual Arrhenius factor in the classical limit when /iAo(rio(Ao) + ■^) ^ fc^T 
for hXo <C ksT This term takes care of the quantum effects of the heat bath which thermalizes the particle inside the 
well. The vacuum term 1/2 prevents the distribution from being singular as one approaches to absolute zero. There- 
fore a significant contribution of quantum correction due to heat bath to the rate enters through both equilibrium 
Wigner function and %. The last term, the quantum depopulation factor is essentially due to the weak coupling of the 
c-number unstable normal mode with the stable modes at the barrier top. This quantity signifies the nonequilibrium 
nature of the dynamics at the barrier top and depends on the average energy loss of the c-number unstable mode as 
well as on the energy dispersion. Since both of them are sensitive to the coupling , we emphasize that the quantum 
effect due to nonlinearity of the system makes its presence felt in these quantities. Apart from these couplings, quan- 
tum nature of the system is also manifested in tp, the round trip time, for a complete traversal in the well. The latter 
point can be understood more explicitly as we go over to Section VI. 

Before closing this section we point out that (i) because of low temperature quantum correction to Kramers'-Grote- 
Hynes factor, the quantum nature of the unstable mode, and a Wigner description of probability distribution inside 
the well, the rate expression is valid both above and below the cross-over temperature. Thus this works well below 
the socalled activated tunneling regime down to vacuum limit, (ii) Since in this present calculation the dissipation 
effects are not restricted to Markovian limit, the rate expression is valid for arbitrary frequency-dependent friction, 
(iii) Finally from the expressions for the two primary quantities that determine the quantum depopulation factor , i.e. 
, the average energy loss AE and energy dispersion {SE'^) it is apparent that the classical limit of the depopulation 
factor can be recovered by reducing these two quantities in the limit fcsT 3> hXo and neglecting quantum correction 
in Vat. The result is the classical depopulation factor 0. (iv) An important advantage of the present c-number 
scheme is that it allows us to use a Gaussian kernel (14.1611 which is exactly classical in form (as used by Melnikov and 
Meshkov) but quantum mechanical in its content since the energy loss AE and the dispersion [{SE'^) = D] are quantum 
mechanical in character. The quantum-classical correspondence thus be immediately restored. The rate expression 
()4.21(l reduces to its classical non-Markovian expression at high temperature where hx{n{x) + 1/2) becomes ksT so 
that the Wolynes factor x goes over to unity and the depopulation factor takes its classical value. We obtain the 
expression for PoUak-Grabert-Hanggi as follows: 



where AE as given by (|4.11l) quantum corrections arising out of system nonlinearity. 

V. THE QUANTUM CORRECTION TO ENERGY LOSS OF THE UNSTABLE MODE 

The key issue in Kramers' turnover problem is the calculation of energy loss AE of the unstable mode during its 
round trip in the well over a time tp. This is extremely sensitive to coupling of the stable and the unstable modes 
described in terms of the equations of motion for the normal modes (|4.4I) and 1)4. 5|l . The present calculation by virtue 
of considering the quantum system mode takes care of this loss through C(i) or through the nonlinearity of the system 
potential V^r. Furthermore the time tp depends explicitly on the quantum corrections as may be seen in the next 
section. Thus the quantum corrections i?„ (in Ea. H2.25|l ) are the key quantities that need to be determined. In an 
illustrative example of the next section with a cubic potential, we show that B2 is the relevant quantity for the leading 




(4.22) 
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order quantum correction. In the following we give a recipe for calculation of quantum corrections. We return to the 
operator equation (|2.2II and put (|2.19|) and H2.20|l to obtain 



5q — Sp 



(5.1) 



sp + j{t- t') Sp{t') dt' + V"{q) 5q 



E i v^^^+'\q) [sr - {sr)) - F{t) - fit) 



(5.2) 



n>2 



We then perform a quantum mechanical averaging over bath states with nt=i{|c^i(0))} to get rid of the term 
F{t) — f{t). The Eas. H5.1|l and H5.2|) along with (|2.1t)|) and H2.17|l form the key element for calculation of the quantum 
mechanical correction due to the nonlinearity of the potential. Considering the friction kernel "f{t) to be arbitrary 
(but decaying) we may calculate the leading order quantum correction for the harmonic mode at the barrier top for 
which higher derivatives of V{q) in H5.2|l vanish. The above equations can then be solved by Laplace transformation 
technique to obtain 



Sq{t) = Sp{0) a(t) + 1 + u;l / C,{t')dt' Sq{0) 



and 



where 



C,it) = L- 



s 7(s) 



and 



Cq{t) = i + ujI a,{t') dt' 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



and 7(s) is the Laplace transform of "f[t) defined as 7(5) = ^*dt. After squaring and quantum mechanical 

averaging Eq. H5.3|l yields 



(Sq'it)) 



)cl(t) 



{5p\0))Cl{t) + {5e{0) 

+ a(i) c,(i)(<5p(o) sm + 6m (Sp(o)) 

The relevant quantum correction in Vn in the leading order is B2 obtained as a time average of {Sq'^{t)). 



(5.7) 



Bo 



{Sf{t)) dt 



P JO 



(5.8) 



While for Markovian friction, the above equation can be calculated analytically , one must take resort to numerical 
evaluation of Cy (t) and Cq (t) for arbitrary friction kernel. We emphasize that the correction B2 is the leading order 
quantum correction for the unstable mode and the friction kernel is considered to be arbitrary in nature. Thus the 
quantum correction due to the unstable mode affects {SE'^) as well as AE through tp and C(^)- The quantum nature 
of the heat bath on the other hand is taken care of through the width parameter of the Wigner canonical thermal 
distribution function. It is pertinent to mention at this point that the quantum correction to the average energy loss 
and dispersion of energy can be calculated by including higher order contribution with the help of the basic equations 
(|5.1|) and (|5.2|) within the framework of the theory. However, the coupling of the unstable and the stable modes at 
the barrier top being weak it is sufhcient to consider this leading order quantum correction for the present treatment. 
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VI. APPLICATION TO A CUBIC OSCILLATOR 

The theory developed so far is fairly general in the sense that it takes into account of an arbitrary form of nonlinear 
potential V{q) and a frequency dependent friction. We now consider a simple nonlinear potential of the form chosen 
for comparison with standard result [Tol | 

V{q) = -aq'-bq^ + ^ (6.1) 

The extrema of the potential corresponds to g = and q = (—go) = ~2b/3a, the respective potentials being 
4b^/27a'^ and respectively. Thus we have 

the metastable minima is at q = —qo- Furthermore we assume a Lorentian form of spectral density function, J{uj) 

as 



J(u}) = Tp-^ — n (6-3) 



for which the friction kernel has an exponential form as 



The Laplace transform of 7(t) is given by 



7(t) = - exp(- — ) (6.4) 

Tc ITc 



7(s) = 7/(1 + s 7 Tc) (6.5) 



We follow the classical procedure of Straub, Borkovec and Berne [T^ll^ to consider the potential ()6.1|l in the form 
of a piecewise continuous harmonic potential as 

V{q) « i (q + q^f for q < ^q* (6.6) 

« Vi-^ Lol q' for q > -q* (6.7) 

From the continuity of the potential V{q) and its derivatives at q = —q*, i.e., 1/2 ojq {q + qo)^|-(j' = Va ~ 
1/2 Lul and Wg (5 + qo)\-q* = -u>l q\-q* we obtain 

= ^ and q*^~- (6.8) 

Now the nonlinearity of the potential around q = can be estimated from the general expression H2.28|l and (|6.1|l . 
Thus we have 

Vn = V2iq) + V3iq) (6.9) 

n>3 n>2 

= -a q^ - 3B2 a q - B3 a (6.10) 

where B2 and B3 defined in (|2.25|l , are the quantum corrections to the potential due to nonlinearity. The details of 
the evaluation of _B„ are given in Sec.V. Linearizing Vn around q = —go we obtain the nonlinear part of the potential. 
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VNiq) ~ {3 a qa) q"^ + {3 a - 3 a B2) q + {a - a S3) 



(6.11) 



and its derivative 



Vpfiq) = {6 a qo) q + {3 a ql - 3 a B2) 



(6.12) 



We are now in a position to write down the equation motion H4.4|l for the unstable mode explicitly for the potential 
concerned for the present problem 



452 

(4 b uIq) p + ( — - 3 a B2 



(6.13) 



or in the following form 



(6.14) 



where Aq = 46 - A^ 



(6.15) 



'46^ 

and r — — [ — — Mqq — 3 a B2 

oft 



(6.16) 



If the unstable mode start moving at — p* with total energy V'^ ^ then the initial value of the momentum for the 
unstable normal mode is 



p{t - 0) 



uoa \ 3a 



The corresponding unstable co-ordinate is 



(6.17) 



p(t = 0) = 

For a round trip, the time elapsed is tp] we have 



- (- 

uoo \3a 



(6.18) 



pit :^tp) = - 

and p{t = tp) 

The equation for the unstable mode H6.13|l can be solved to obtain 



-(- 

uqo \ 3a 



(6.19) 
(6.20) 



pit) = M cos(Aot) + N sin(Aot) + 



where M and N are given by 



(6.21) 



3uo() a Aq / AoMoo V ■^^ 



(6.22) 
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The round trip time tp can be calculated form Eq. (|6.21(l and its derivative equation for p and applying the conditions 
(|6.19|) and H().20|l on them so that we have 

AP - iV^ 2 M N 

cosiXotp) - J^t:^ , sin(Aotp) = j^TJ^ (6-23) 

Here we also note the range of tp as n < Xatp < 2tt. An important point to emphasize is that tp contains a quantum 
correction. 

In order to calculate the energy loss due to the unstable mode we now return to Ea. ()4.11|l so that we write it as 

A£; = i [ " dt [ " dt' Kc{t - t') C{t) C{t') (6.24) 
2 Jo Jo 

where Kc{t) has the form given by Ea. H2.46|l since = (wio/^oo)^ 

Since 7(5) is given by (|6.5|l . Kc{t) can be calculated explicitly using (|2.47|l as in the classical theory to obtain 



K,{t) = i exp(-ei) (1 + 2e) cosh((7t) + + 2e)^ h ^-^^^^^-^ 

exp{-Xbt) (6.25) 



2 

1 
2 



where ^ and a are defined as ^ = ^ (Af, + l/{jTcj) and tr^ = — Lul/{jTcXb)- C{t) is given by (|4.6fl and using (|6.12fl 
we obtain explicitly 



'46^ 

at) = -4 b ul. Pit) - i-^ulo - 3aB2 j (6.26) 
We then make use of the solution 16.21|l for p{t) in (|6.26|) to obtain 

C(t) = Ml cos(Aot) + ^1 sin(Aot) - P (6.27) 

where 

Ml = -4M6uoo 

Ni = -4:Nbulo (6.28) 

4 6 u^n r 4b'^u?,n , 
■ — 3ah2 



^ - -^00 ' u-oo „ „ 



We are now in a position to calculate the energy loss AE of the unstable mode using (|6.24() from (|6.25|) and H6.27|l . 
After a little bit of algebra we obtain 

AE = [i?(e + <j) + i?(e - <j)] 

- + ^'^j - [fl(g + ^) - - a)] - 1 E(A,) (6.29) 

where 

i?(2) = ["dtf^ di'e~^(*~*') [Ml cos(Aot) + A^i sin(Aot) - P] 
Jo Jo 

X [Ml cos(Aot') + ^1 sin(Aoi') - ^] (6-30) 

Explicit evolution of the integrals in H6.30|l yields 
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Riz) 



P Miz - NiXo 

z 



z [ Xl + z^ 



cos(Aoip) 



(MiPXq + MiNiz - PNiz Xoz{Ml + TVf) 

{Xl + z^Y 
2P{Miz - NiXo) 



V 4^1 + ^') 

T"^ 2{Xl + z^) . 



sin(Aotp) ^ e 



Mfz^ _ ^2^2 7Viz(Mi - 4P) P2 



(6.31) 



Our next task is to calculate the width D of the distribution function l|4.16|l . D is defined in H4.14|l (and (|4.13|) 'l 
and we now rewrite {SE'^){= D) as 



where 



(,5p2) = r dt r dt' KQ{t- 1') at) at') 

Jo Jo 



(6.32) 



We also recall Eq. 



N 

^^Q{i) = '^af ( + § ) cos(A,t) 

2—1 



(6.33) 



N 



Kcit)^Y.9'i cos(A,t) 



i=l 



By going over to continuum modes Kc{t) and KQ{t) can be expressed as 



and 



KQ{t)^ / dAg^(A)/(A)J(A)cos(At) 
Jo 



K^{t)^ dXg^{X)J{X)cos{Xt) 



(6.34) 



(6.35) 



(6.36) 



where /(A) = ^ fi,A coth(ftA/2fcBT). It can be easily shown that the above two quantities can be related in the 
Fourier transform domain as 



KgiX) - f{X)K,iX) 

Here Kq{X) and Kc{X) are the cosine transform of KQ{t) and Kc{t) respectively. 



(6.37) 



KcW = - 



1/2 



dt Kelt) cos(Ai) 



(6.38) 



Making use of the expression Kc{t) from Eq. 
The inverse transform of Kq{X) finally gives 



Sj) in Eq. 



followed by a multiplication of /(A), yields Kq{X). 
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KQit) = 



Here 



Wi 
W2 



1/2 



Wi 



dX -?iAcoth 
2 



hX 



cos(At) 



2kBT 
W2 W3 



(^ + a)2 + A2 (^_a)2 + A2 A^ + A^ 



2^ 


r 










2^ 


r 










2^ 


r 


Xb 


71" y 




2 



l + 2e (l + 2e)e-A6 



4ct 



4(7 



with H6.27(l and Hfi.39|l we obtain the expression for {6E'^) as follows. 



(6.39) 



(e + 'T) 



(6.40) 



(5£;2) = L> 



1 



1/2 



n / dA Acoth 
2vr/ Jo \2kBT ^ 



hX 



(C + a)2 + A2 (e-fT)2 + A2 A2 + A2 

X [" dt [ " dt' cos[A(i - i')]C(OC(^') 
Jo Jo 



(6.41) 



The quantum expression of the depopulation factor in Eq. H4.21|) primarily depends on the average quantum energy 
loss AE and dispersion (SE^). The expressions of these two quantities have been given in Ea. H6.29|) and Ea. (|6.41|l for 
a cubic oscillator (|6.1|) . An explicit evaluation of these quantities for calculation of the rate or depopulation factor 
requires numerical calculation. The parameter space chosen for illustration of the main results depicted in Figs. 2-5 
is based on the earlier work by PoUak et al Rips and Pollak [13, Straub et al |l3l | for a comparative study. The 
results are given below. 

In Fig. 2 we exhibit the quantum turnover by plotting the quantum rate, F/Fg, where Fq = 



(ti;o/27r)xexp — {n\\ )+!) ' ^ ^ function of dissipation parameter 7 at three different temperatures kgT = 5.0 

(dotted line) and ksT = 1.0 (dashed line) and ksT — 0.0 (solid line) for a model potential with = 20.0, 
a = 0.03042, b = 0.5. The calculation here been carried out with a non-Markovian friction kernel for Tc = 1.333 
and the quantum correction has been taken upto second order. In order to allow ourselves a fair comparison with 
classical non-Markovian theory of Pollak, Grabert and Hanggi, we compare the quantum (solid line) curves with the 
corresponding classical (dotted line) curves at two different temperatures ksT = 1.0 and fcsT = 2.0 in Fig. 3 for the 
same set of parameter values as in Fi g. 2. The following two points are noteworthy . We observe (Fig. 2) in confor- 
mity with the earlier observation |llL l23l| that as the temperature is lowered the maximum at which the quantum 
turnover occurs shifts to the left and the damping region that corresponds to classical energy diffusion regime becomes 
exponentially small as one approaches to absolute zero. Second, we observe that the strong friction tends to make 
the dynamics more classical since the quantum correction is suppressed by dissipation in this regime. On the other 
hand differential behaviour of the rate in the classical and quantum regime is felt at weak friction regime. This is in 
confirmation with earlier observations. In Figs. 4 and 5 we show a comparison of the quantum rate, F/Fg, calculated 
on the basis of the present theory (solid line) with that of Rips and Pollak (dotted line) for — 26.188, a — 0.0266 
and b — 0.5 at two different scaled temperatures fcsT = 1.788 and fc^T = 2.617 respectively in the Markovian limit. 
The agreement is found to be quite satisfactory. Thus although the depopulation factor in the rate expression (|4.21|) 
apparently differs in form from that of Pollak and Rips based on quantum transition probability obtained from the 
solution of master equation, the good numerical agreement between them can not be overlooked (since the remaining 
factor in F/Fg, the Grote-Hynes factor is same in both the expressions). Figs. 4 and 5 therefore serve as a consistency 
check of the present calculation. 
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VII. SUMMARY AND CONCLUSION 

Based on a quantum Langevin equation we have constructed a c-number Hamiltonian for a system plus N-oscillator 
bath model. This allows us to formulate a normal mode analysis to realize a c-mimbcr version of the multidimensional 
transition state theory and to derive a quantum expression for the total decay rate of metastable state. The result 
is valid for arbitrary damping strength and noise correlation and temperature down to vacuum limit. The theory is 
illustrated on a cubic potential with non-Ohmic dissipation. The following pertinent points are noteworthy. 

(i) We have shown that the expression for quantum rate coefficient is a product of five terms , e.g. , classical well 
frequency, Kramers'-Grote-Hynes factor, a vacuum corrected or generalized Wolynes factor representing quantum 
transmission and reflection, an exponential term corresponding to Wigner canonical thermal distribution i.e. the 
generalized Arrhenius term and a quantum depopulation factor. Of these the Wigner term and the vacuum corrected 
Wolynes term refer to equilibriation of quantum particle in the well and therefore corresponds to a c-nunibcr multi- 
dimensional transition state result. Since the distribution unlike the Boltzmann is valid even as T ^ 0, the quantum 
eff'ect due to heat bath can be well accounted by these terms even below the activated tunneling regime. 

(ii) The quantum depopulation factor has a form which is very much similar to its classical counterpart, although 
intrinsically the relevant quantities determining this factor, i.e. , the average quantum energy loss and dispersion are 
quantum mechanical in their content. 

(iii) The classical limit of the quantum rate expression depends on Wigner, generalized Wolynes and depopulation 
factors. It is easy to see that they reduce to Arrhenius factor, unity and classical depopulation factor, respectively in 
the limit fku <^ ksT . 

(iv) Since the quantization of the nonlinear system mode adds an additional contribution to depopulation factor 
arising out of the weak coupling between the unstable mode and the stable modes the rate is significantly modified. 

(v) The theory takes into account of the quantum eiT(X"ts in full due to heat bath, while the quantum correction due 
to nonlinearity of the system potential can be calculated systematically order by order to a good degree of accuracy. 
In the illustrative example with cubic potential we have considered the quantum correction of the second order. 

(vi) The quantum depopulation factor interpolates the energy diffusion limited regime to spatial diffusion regime 
and depicts the correct turnover scenario down to absolute zero. The maximum of the turnover shifts to the low 
damping regime as one approaches the absolute zero. The quantum rate is significantly enhanced over the classical 
rate in the weak damping regime while strong dissipation overwhelms the quantum correction. 

(vii) An important aspect of the present theory is that it takes care of activation and tunneling within a unified 
description and is equipped to deal with the rate at temperature down to vacuum limit. This is a distinct advantage 
over path integral Monte Carlo method since numerically the relevant propagator poses serious problem as the 
temperature approaches absolute zero. 

(viii) The present theory is a synthesis of the classical formalism of normal mode analysis and weak coupling theory 
of unstable and stable modes for calculation of average energy loss and dispersion within the framework of our recently 
developed c-number quantum Langevin equation. It is important to emphasize a distinct advantage of the present 
treatment. The present c- number formalism allows us to formulate the Guassian transition probability that the energy 
of the unstable mode changes for one period of motion within a " classical- like" prescription, rather than by solving a 
quantum master equation restricted to a Markovian description. 

The quantum turnover theory as presented here is based on a canonical quantization procedure and positive definite 
Wigner's thermal distribution for harmonic oscillators rather than path integral or master equation formalism. It can 
be readily applied to other models t e.g. , a double well oscillator and to the problems of dephasing and related issues. 
The systematic improvement can be made by taking care of the quantum corrections of higher orders. 
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Figure Captions 

Fig.l. A schematic plot of the potential defined in Eq.(6.1) 

Fig. 2. Quantum rate, F/Fq, (where Tq — {llJo/2tt)x exp — ), in the non-Markovian regime is plotted 

against dissipation parameter, 7, for three different temperatures fc^T = 0.0 (solid line), ksT = 1.0 (dashed line) and 
ksT = 5.0 (dotted line) with = 20.0, a = 0.03042, b = 0.5 and = 1.333. 

Fig. 3. Quantum turnover (solid line) is compared with classical turnover of PoUak, Grabert and Hanggi (dotted 
line) in the non-Markovian regime for two different temperatures ksT — 1.0 and fc^T = 2.0, other parameters remain 
same as Fig. 2. 

Fig.4. The comparison of Kramers' turnover for = 26.188, a = 0.0266, b = 0.5 at kT = 1.788. The solid Une 
represents the quantum turnover by the present method, the dotted line represents the calculation by Rips and PoUak 

m 

Fig.5. The same as in Fig.4 but for ksT = 2.617 
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